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AC conductivity spectra of xAgI-(l-x)AgP03 fast ion conducting glasses spanning the flexible, 
intermediate (isostatically rigid) and stressed rigid phases are analyzed. The rescaled frequency 
dependent spectra are mapped into time-dependent mean square displacements out of which a 
typical lengthscale characterizing the spatial extent ^/(i? 2 (oo)} of non-random diffusion paths is 
computed. The latter quantity is studied as a function of Agf composition, it is found to display a 
maximum in the intermediate phase, providing the first clear evidence of a typical lengthscale of a 
dynamical nature when a system becomes isostatically rigid and enters the intermediate phase. 
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An amorphous network progressively stiffens and be- 
comes rigid as its connectivity or mean coordination 
number f increases. From a mechanical viewpoint, such 
an evolution can be understood using rigidity theory 
which considers nearest-neighbour interactions acting at 
the microscopic level |l|, and enumerates the average 
number of constraints n c per atom and its balance with 
respect to the atomic degrees of freedom (3 in three di- 
mensions). This has led to the recognition of a rigidity 
transition [2|, |3j separating underconstrained (or flexible 
when n c < 3) from overconstrained networks (stressed 
rigid, when n c > 3). Numerous experiments have con- 
firmed these simple predictions 4] , especially in glass sci- 
ence where bulk chalcogenide and oxide glasses have been 
studied in detail. 

While the original theory predicted a single optimized 
glass composition where the microscopic structure is iso- 
static (having n c ~ 3), more recent experiments have re- 
vealed a second transition [5] for various network-forming 
glasses, providing a finite width, offering now a whole 
range of such isostatic compositions, and defining an in- 
termediate phase (IP) bounded by the flexible (at low f 
or n c ) and stressed rigid phase (at high f). The signature 
of this phase has been detected in optical [8| , calorimet- 
ric [6j, and electrical [7( with clear changes in behaviour 
when the system becomes isostatically rigid for e.g. com- 
plex heat flows at the glass transition or Raman optical 
elastic power-laws. Yet no typical spatial quantity giving 
rise to a characteristic distance (or lengthscale) has been 
observed in the IP, although such an isostatic lengthscale 
has been predicted to grow [•>] in jammed soft sphere 
solids which bear some similarities to network glasses. 
Neutron and high energy X-ray diffraction on chalco- 
genide glasses have been recently reported [lfj, 111], al- 
though a structural origin of the IP in typical quantities 
(e.g. position, width and height of the first sharp diffrac- 
tion peak) related to the static structure factors was not 



observed. 

In this Letter, we show that a lengthscale associated 
with the intermediate phase appears in the fast ionic 
conductor (l-x)AgPC>3-xAgI which has an intermediate 
phase in the 9% < x < 37% Agl range [l2| . The present 
quantity characterizes the spatial extent of subdiffusive 
ion-motions, and is computed from frequency-dependent 
conductivity and permittivity data together with linear 
response theory. The lengthscale appears to be strongly 
correlated with calorimetric, optical or electrical trans- 
port quantities that usually reveal the IP, but is nei- 
ther associated directly with static network structure nor 
with changes in molar volume, as both evolve smoothly 
with glass composition. Another lengthscale character- 
ising the typical distance mobile ions travel to overcome 
backward-forward driving forces is found to be only sen- 
sitive to the onset of network flexibility at x=37% Agl 
[121 ]. The present findings highlight, therefore, the clear 
dynamical nature of the lengthscale underlying the IP, 
and underscores a possible connection with the peculiar 
relaxational phenomena of isostatic glass-forming liquids 
0. 

Details on the synthesis and the structural, calorimet- 
ric, electrical and spectroscopic properties of dry and ho- 
mogeneous bulk (l-x)AgP03-xAgI glasses were discussed 
earlier in Refs. [14||-[15l|. Conductivity measurements 
were performed on disks about 10 mm in diameter and 
2 mm thick on which Pt electrodes were deposited. The 
complex impedance Z*(w)=Z'(u;)-r-iZ (w) was measured 
by a Solartron SI 1260 impedance analyzer over the fre- 
quency range of 1 Hz-1 MHz. Three temperatures have 
been considered: 250 K, 300 K and 350 K. Conductivity 
a*(uj) was deduced from the complex impedance Z*(w), 
sample thickness, t, and surface area S covered by plat- 
inum using the formula: cr*=t/Sxl/Z*. The real part of 
the dielectric permittivity e'(u) was computed using the 
electrical modulus M*: M*^l/e*=icueo/a* , where e is 
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FIG. 1: (color online) (a) Experimental conductivities of a 
25AgI-75AgP03 glass at three different temperatures and cor- 
responding (c) dielectric permittivities. The red horizontal 
line serves to define the low frequency permittivity e'(0) used 
in equ. p)l. (b) Master curve a fade as a function of the 
rescaled frequency v fodcT for three selected compositions at 
T=250 K, derived from the experimental conductivities. 



the permittivity of free space. Because of electrode po- 
larisation occuring at higher temperatures (see Fig. lc 
and text below), we have restricted the analysis to the 
temperature of 250 K. 

Figure QJ, shows the log-log plot of the real part of 
the conductivity with frequency at three different tem- 
peratures for a glass at x=25% Agl. Dielectric permittiv- 
ity spectra e'(w) were computed (Fig. [TJ;) and tracked 
with composition (jj. Concerning the latter, electrode 
polarization effect onsets at low frequency manifested by 
a rapid growth of e'. This phenomenon hides the in- 
trinsic frequency behaviour of s'(u>) so that an evalua- 
tion of the limit e'(0) (defined in Fig lc) is only pos- 
sible at low temperature [lfjj. Similar spectra and fea- 
tures were obtained at other compositions and tempera- 
tures (see Ref.[7]). At low frequency and low temper- 
ature (250 K), we obtain a dc regime which displays 



an Arrhenius behaviour, whereas at higher frequencies 
a dispersive regime is observed, with a or e' becoming 
frequency dependent. A frequency i/ p (x,T) (also called 
dielectric loss peak frequency) satisfying a{v p )/adc=2 is 
usually introduced to characterize the cross-over range 
between the dc and dispersive conductivity regimes, the 
latter being associated with a subdiffusive regime in the 
time domain valid for t<tp=l/(27ri/p). The frequency 
Up increases with increasing temperature, and a simple 
rcscaling by a factor ad c T has been proposed |l2j in or- 
der to produce a master-curve for various temperatures 
and compositions. Results are shown in Fig. lb for three 
selected compositions which show indeed that the data 
nearly map onto the same curve. 

We now build on the linear response theory developed 
by Roling and co-workers I la , Il9j and recently applied 
to borophosphate glasses [20j. It allows extracting from 
conductivity spectra a mean square displacement (r 2 (t)) 
of the mobile ions given by: 



(r 2 (t)) = 



12k B TH R 
N v qH 
= (R 2 (t))H R 



dt' 



a{v) 



sin(2irvt')di> 



(1) 



N v is the charge density, and (R 2 (t)} is the mean-square 
displacement of the center of charge of the mobile ions. 
H R is the so-called Haven ratio [191 ] . found usually be- 
tween 0.2 and 1.0 and depends on the concentration of 
charge carriers 2l|, |22j. It characterizes the degree of 
cooperatvity of ion motion and can only be obtained by 
combining measurement of <7d c with the low-frequency 
limit of the self-diffusion constant -D'(O). As Hr is not 
known for the present system, we focus on ( R 2 ( i)). From 



the permittivity spectra, one can extract [18lll9j | the long- 
time limit of a rescaled quantity, (i? 2 (oo)}, using (R 2 (t)), 
given by: 



(i? 2 (oo)> 



limt-t-oo-^i 



6k B Te 



N v q 2 



(R 2 (t))-6D'(0)t 

e'(0)-£'(oo) 



(2) 



Figure 2 shows (R 2 (t)) for three selected compositions 
corresponding to the flexible (40% Agl), intermediate 
(25%) and stressed rigid phase (4%), computed from 
equ.(l). At long time (t>t p ), (R 2 (t)) displays the diffu- 
sive regime as detected from the slope of 1 in the log-log 
plot, while a sub-diffusive regime corresponding to corre- 
lated forward-backward motions appears at shorter times 
scales (t<t p ). Here the limit t p (or its rescaled quan- 
tity tpUdcT) designates the time beyond which ions start 
to diffuse. One can therefore consider the characteristic 
length ^J{R 2 {t p )) to be the typical distance mobile ions 
have to travel to overcome the interactions responsible for 
these forward-backward motions (see Fig. [5] inset). Ac- 
cording to the definition of the self-diffusion coefficient, 
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FIG. 2: (color online) Mean-square displacement (R 2 (t)} at 
250 K of (l-x)AgP03-xAgI at three selected compositions 
representing the three phases of interest, derived from the 
master curve represented in Fig. 1 and equ. (1) as a func- 
tion of rescaled time tadcT- The vertical broken line corre- 
spond to the abcissa t p ad c T defining t p for a glass at x=4%, 
the approximate limit between the diffusive and sub-diffusive 
regimes. Inset: Trajectory of an ion between an initial and 

a final state, which serves to define the quantities \J (R 2 (oo)} 

and y/(B?(t p )). 



for mobile ions the quantity (r 2 (t)) — 6D'(0)t vanishes at 
the (diffusive) long-time limit. This property is not sat- 
isfied for the related (R 2 (t)) characterizing the motion of 
the center of charge of the mobile ions, and the quan- 
tity (R 2 (oo)) given by equ. (2) provides a measure of 
the spatial extent (Fig. 2, inset) of sub-diffusive motions 

M. 



Both lengthscales \/(R 2 (t p )) and J (R 2 (oo)) can be 
determined from the conductivity and permittivity spec- 
tra, and we follow these with Agl content of glasses. Re- 
sults are displayed in Fig. 3. 



(R 2 (co)) displays a max- 
imum in the IP, driven mainly by a maximum in the per- 
mittivity difference e'(0)— e'(oo) appearing in equ. (2) Q. 
Starting from 5.5 A at x=0, this lengthscale increases up 
to a maximum of nearly 9 A in the IP and then decreases 



to 2.5 A in the flexible phase. Trends in \l (R?(oo)) as 
a function of x clearly correlate with those in the non- 
reversing heat flow, AH nr , obtained from calorimetric 
measurements [121 ] . and with the three regimes for dc 
conductivity. These trends cannot be explained on the 
basis of free volume changes induced by Agl doping as 
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FIG. 3: (color online) (a) Plot of the characteristic length- 
scales a/ (R 2 (oo)) (blue) and \/(R 2 (t P )} (red, right axis) in 
(1-x) AgP03-xAgI as a function of Agl composition, compared 
to (b) (extracted from Ref. [7]]): conductivity aa c (black 
filled symbols) at different temperatures and non-reversing 
heat flow AH nr (open red circles, right axis). The vertical 
broken lines define the intermediate phase boundaries. 



molar volumes are known [15| to decrease linearly. These 
results suggest that there is a maximum in spatial extent 
of subdiffusive motions in the IP. At low Agl concentra- 
tions, isolated subdiffusive regions (SDR) are embedded 
in a stressed rigid structure and lead to a low spatial 
extent. The increase of the number of hopping sites in 
the IP 7] leads to a growth of the subdiffusive regions 

(and y(R 2 (oo)}). Finally, once these SDR percolate, 
as also suggested from Molecular Dynamics simulations 
[23j, the spatial extent of these regions decreases. Fig. 
3 illustrates for the first time that a typical lengthscale 
is associated with the IP, which is of dynamical origin. 
There is no structural origin of such a lengthscale be- 
cause relevant quantities deduced from X-ray diffraction 
24] behave smoothly with composition, an observation 



11] 



also reported on a different systems in 

The other quantity of interest, the lengthscale 
J (R 2 (t p )) , is found to display a threshold at the inter- 
mediate to flexible transition (Fig. 3, right axis), also 
known as the rigidity transition. In fact, as long as the 



system is rigid (stressed or isostatic) , this lengthscale re- 
mains nearly constant at about 0.5 A, largely because a 
large energy is needed to locally deform the network in or- 
der to overcome the backward driving forces causing the 
correlated subdiffusive forward-backward motions. One 
control parameter of y/(R 2 (t p )), the time t p for the on- 
set of diffusion decreases from 400 /is at x=0 to 10-15 /is 
close to the flexible transition (37-40%). Once the sys- 
tem has become flexible and floppy mode proliferate, ion 
motion is facilitated and the time decreases to t p =6.7 /is 
at x=45% Agl. At this concentration, and because of the 
large number of carriers, the time needed to leave a SDR 
reduces, but a larger distance is needed to overcome the 
backward driving forces induced by the presence of an 
increased number of Ag cations. 

In summary, we have shown for the first time that a 
typical lengthscale, of a dynamical origin, is associated 
with the Intermediate Phase. This lengthscale appears 
in fast-ion conducting glasses, and measures the spatial 
extent of subdiffusive ionic motions that increase when 
the system becomes isostatic and decreases once conduc- 
tion pathways percolate in the flexible phase. Another 
lengthscale characterizing the forward-backward motion 
at play in the subdiffusive regime is found to signal only 
the intermediate to flexible transition when the network 
softens and leads to a strong increase of the dc conduc- 
tivity. These results found for a fast ion conducting glass 
may also exist in covalent network glasses. Anomalies 
in viscosity have been observed in network glass-forming 
liquids such as binary Ge-Se [25j , and since viscosity and 
diffusion are related via the classical Nernst-Einstein re- 
lation, such a dynamical lengthscale may exist in e.g. Ge 
and Se diffusion measurements. 
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